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Abstract: The magic couple of SUSY and GUT still appears the most elegant and predictive 
physics concept beyond the Standard Model. Since up to now LHC found no evidence for super- 
symmetric particles it becomes of particular relevance to determine an upper bound of the energy 
scale they have to show up. In particular, we have analyzed a generic SUSY-GUT model assuming 
one step unification like in SU(5), and adopting naturalness principles, we have obtained general 
bounds on the mass spectrum of SUSY particles. We claim that if a SUSY gauge coupling unifica¬ 
tion takes place, the lightest gluino or Higgsino cannot have a mass larger than ^ 20 TeV. Such a 
limit is of interest for planning new accelerator machines. 
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1 Introduction 

The recent detection in July 2012, of a new particle compatible with the Higgs boson by ATLAS 
[1] and GMS [2] collaborations at LHG certainly represents a milestone that has once again con¬ 
firmed the predictive strength of the Standard Model (SM) of strong and electroweak interactions. 
However, from the theoretical point of view it is hardly to believe that the SM is the last step 
toward the unification in a simple principle of all the fundamental interactions. Several phenomena 
or open problems suggest the presence of physics beyond the SM. One can just remind few of them, 
like the hint of unification of all gauge couplings for extreme large energy, the particular structure 
of fermion masses, the baryonic/leptonic number violation processes necessary at scale larger than 
the electroweak one in order to yield baryogenesis, the problem related to the separation of very 
different energy scales in a field theory with scalars {hierarchy problem), etc. 

The Grand Unified Theory (GUT) paradigm is able to address part of these problems. Sug¬ 
gested by the simple extension of the SM scheme to larger gauge groups, it has represented a huge 
scientific effort, which, starting from seventies, has had to face with the experimental counterpart 
of proton decay search, and with the necessity of an even too rich scalar content of the theory. The 
main feature that makes plausible a GUT scheme at very high energy, namely a compact gauge 
group describing the fundamental interactions, comes from the observation that running gauge 
couplings tend to get closer and closer with the increase of energy. 

A GUT is generally characterized by a typical energy scale, Mqut- It is naively defined as 
the scale where the SM gauge couplings cross at one point according to their running. The same 
definition can be rephrased as follows: the GUT scale is the scale at which a compact gauge 
symmetry G spontaneously breaks down to SM (namely SU{3) x SU{2) x U(l)). Below such a 
scale the gauge couplings, being equal to each other at /r = Mgut, run down and match their 
experimental values at the benchmark point p, = Mz- However, this is strictly true for SU{5) 
model only, for which both the above definitions are equivalent. Moreover, in 517(5) the great 
desert is occurring since there are no intermediate thresholds due to new massive states starting 
from Mt up to GUT scale. 

Then a question is in turn and concerns the size of Mqut- While a natural upper bound is 
represented by string or reduced Planck scale Mp ~ 2 • 10^® GeV, the lower limit is settled by the 
predicted degree of stability of ordinary matter (proton decay). The experimental limits on the 
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proton decay mediated by the heavy gauge bosons pose stringent bounds on their masses, hence on 
Mgut- 

Before LEP epoch, non-supersymmetric GUTs could not been discriminated due to the poor 
precision characterizing the measurements of SM parameters, even though the low energy scale of 
unification, order < 10^® GeV, was already in strong tension with the proton lifetime limit. Such 
situation then drastically changed with the precise measurements at LEP, which essentially excluded 
non-supersymmetric GUT. In fact, after the precise measurements of the gauge couplings Oi, 02 , 
03 at Mz, and of other SM quantities at LEP, it was clear that a real one step unification of running 
gauge couplings was not possible at high energy, at least without assuming new physics. However, in 
the framework of Supersymmetry (SUSY) one obtains a milder running of couplings that provides a 
satisfactory unification at a common high energy scale Mgut [3-7]. Moreover, SUSY also provides 
a road map for solving the gauge hierarchy and Doublet-Triplet splitting problems. In particular, 
minimal SU(5) allows a technical solution for the hierarchy and Doublet-Triplet splitting via a fine 
tuning, which is stable against radiative corrections due to SUSY [ 8 , 9]. However, there exist more 
natural solutions without fine tuning, as ‘Missing Doublet Mechanism" (MDM) in extended SU(5) 
[10-14], “Missing v.e.v. Mechanism" (MVM) in SO(IO) [15-18], or “pseudo-Goldstone boson instead 
of Fine Tuning" (GIFT) Mechanism in SU( 6 ) [19-24].^ 

The golden goal of supersymmetry was to explain naturally the origin of the electroweak (EW) 
scale in terms of the soft supersymmetry breaking (SSB) parameters, without artificial fine tunings 
among them. In this design of a natural SUSY, the EW scale has to be originated by the SUSY 
breaking scale itself. This means that already above few hundred GeV the SM has to be replaced by 
the supersymmetric extensions of the Standard Model (MSSM in the minimal version), which would 
imply a rich new phenomenology that LHG should have already found. With the conclusion of the 
8 TeV LHG run I, MSSM have been greatly constrained by a variety of direct searches [34-36] . The 
properties of the discovered Higgs boson are fully compatible with the Standard Model Higgs, and 
a second light Higgs predicted in natural SUSY was not found. In the meanwhile the experimental 
lower bounds on the masses of gluinos and squarks were increased up to TeV scale. 

More in general, the missing evidence up to now for a SUSY phenomenology at LHG has caused 
a pessimistic attitude of the scientific community toward SUSY and SUSY-GUT paradigms, thus 
forcing theoretical physicists to look for different approaches. As a radical possibility, there emerged 
a concept of split supersymmetry [37, 38], which pushes some SUSY partners (squarks and sleptons) 
to arbitrarily large scale while leaving the others (gauginos and Higgsinos) at the low scale. 

In view of this, one can ask if the experimental situation after LHG is really so dramatic to 
justify almost a total give up of the SUSY scheme, or SUSY still remains the most promising and well 
defined concept that particle physics above SM should tend to. In any case one cannot exclude the 
possibility that SUSY lives just around the corner, and gluinos or squarks will be indeed discovered 
at the TeV scale in the second run of the LHG. However, it is important to remark that even before 
LHG, already in the LEP epoch the experimental data supplied mounting evidences disfavoring, at 
some level, SUSY at few hundred GeV [39] . It is worth reminding the problem of electric dipole 
moments of electron and neutron, which within ~100 GeV range for SUSY were predicted too large, 
while increasing the SUSY scale to few TeV they are naturally suppressed to the level of the present 
experimental limits. The situation is similar for the flavor violating processes like /r —> eq, 6 —>■ sq 

^The rudimentary idea of Higgs emerging as pseudo-Goldstone boson was first introduced at the level of SUSY- 
SU{5) model in Refs. [25—27], assuming that Higgs superpotential has an ad hoc global SU{6) symmetry at the price 
of introducing an extra singlet supermultiplet. These models are not fully realistic without local SU{6) completion. 
At the level of local SU{5) symmetry the assumption of Higgs superpotential having SU{6) global symmetry is 
equivalent to many fine tunings, but at the same time such extended symmetry is explicitly violated in the Yukawa 
sector. Let us remark that the supersymmetric pseudo-Goldstone mechanism was also discussed in the context of 
SO(IO) [28], SO{n) and Eg [21]. More recently such mechanism was also applied for the solution of the little hierarchy 
problem in the framework of minimal low scale extensions of the MSSM [29-33]. 
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etc. For arbitrary soft parameters, the flavor violation limits would require the increasing of the 
SUSY scale up to 100 TeV, but such a limit can be lowered to TeV scale by a particular use of 
symmetry arguments, namely by means of an approximate alignment of the soft SUSY parameters 
with the Yukawa couplings. In this concern, one can mention particular theoretical scenarios based 
on a gauged family symmetry SU{3)h [40-43], which can be implemented for SUSY-GUT as well 
[40, 41, 44[. In this case, the flavor violation limits would allow SUSY scales as small as 1 TeV. Later 
on such a possibility was described as a paradigm thereafter simply denoted as Minimal Flavour 
Violation [45[. 

Finally, as will be clarified in the following, the Higgs physics, modulo specific conspiracies, 
seems to indicate that the SUSY scale should be above few TeV as well. Therefore, one can 
conclude that SUSY scale larger than few TeV is quite reasonable, and if SUSY indeed lives near 
this lower bound, then it should show up in the LHC run H. 

In this paper, we try to reanalyze the issue showing the room still remaining for SUSY-GUT 
inspired models under some natural assumptions, which will be clarified in the following, and the 
consequences of such schemes on LHC run H and on future colliders like FCC-ee and FCC-hh 
[46, 47[. Under such assumptions, we study the implications on SUSY mass spectrum due to the 
gauge couplings and Yukawa couplings unification. We adopt the paradigm of one step gauge 
unification, which we simply denote as SU(5) bottleneck, and that still remains a quite general 
condition. GUTs may be based on a group larger than SU{5) like 50(10), SU{6) or Eq, which 
provide promising and predictive models describing physics below the string scale or reduced Planck 
scale Mp. Such latter symmetry groups do not provide automatically one step unification of the 
gauge couplings, and generically they could be broken down to the SM passing through different 
intermediate stages of the gauge symmetry breaking. In this case the gauge couplings unification 
phenomenon would not be a clean prediction of the scheme, but rather an ad hoc reconstructed 
phenomenon. In this respect, the condition of SU(5) bottleneck simply requires that the gauge 
symmetry is broken in such a way that below the scale Mgut it reduces to 5C/(3) x SU{2) x U{1) 
with three gauge couplings ai, 021 0:3 that run down in energy starting from the same value cigut 
( denoting the GUT gauge coupling at Mgut)- In other words we require that independently of 
whatever is the “mother" unified gauge group, and the intermediate breaking scales it may require, 
at the last symmetry breaking stage, nearby Mgut, it behaves as 51/(5). Below such a scale the 
theory should reduce to the Minimal Supersymmetric Standard Model (MSSM). 

In the present analysis we require that all coupling constants, including Yukawa couplings, 
must be order one in the “mother" GUT scenario above Mgut, assuming that there are no ad hoc 
small parameters and no ad hoc fine tuning among parameters. Moreover, the mass parameters 
involved in the SUSY breaking are supposed to be of the same order of magnitude (modulo possible 
differences between F- and Zl-terms). Small Yukawa couplings for light families in the MSSM 
superpotential below the GUT scale do not contradict the previous requirements as we will discuss 
in the next section. 

In summary, to perform our study we require: 

• Unification of all gauge couplings at a single energy scale (Mgut) without intermediate sym¬ 
metry scales {SU(5) bottleneck). We take into account only SUSY particle thresholds, simply 
denoted as SUSY thresholds, at which SUSY particles show up, and GUT thresholds related 
to GUT multiplet fragments that can be lighter than Mgut- Since we do not allow the pres¬ 
ence of ad hoc small parameters these fragments cannot be much lighter than one order of 
magnitude of Mgut- 

• Gonsistency of third family fermion masses: correct mass of the top quark and Yukawa b-r 
unification. 


- 3 - 


• Consistency with the experimental limit on proton decay. 

• The absence of special fine tunings among the parameters (couplings 0(1) at Mgut), which 
can be seen as a general naturalness requirement. 

Moreover, we do not consider the possible limits on SUSY-GUTs coming from the assumption that 
one of the SUSY particle must necessarily be a DM candidate. This is due to the possible presence 
of R-violating terms that would make unstable the lightest SUSY particle at cosmological time- 
scale. Hence, in order to get model independent predictions we prefer do not use such cosmological 
constraints. 

Under this ansatz, we analyze the limits on the mass spectrum of SUSY particles once the 
previous requirements are fulfilled, and indeed we obtain an upper limit of about few TeV in the 
more natural case when all supersymmetric particles have masses of the same order of magnitude. 
The situation changes if one takes into account a possible spread among the superpartner masses 
of different types (gauginos, sfermions and Higgsino). In this case, we look for the minimum of the 
mass spectrum, identifying such value as the energy scale upon which SUSY phenomenology has 
to be detected. By spanning on all compatible SUSY-GUT models (see previous requirements), 
in presence of both SUSY and GUT thresholds, we find that the above minimum can be as large 
as ~ 20 TeV. This represents the main result of our analysis since it provides an almost model 
independent upper bound for the appearance of SUSY phenomenology. 

In order to get such a result, once that all masses of the model are fixed, by using a Mathematica 
code that solves the set of Renormalization Group Equations (RGEs) at two loops, we determine the 
unification GUT scale Mqut and the possible compatibility of such a choice with the experimental 
measurements of gauge couplings at EW scale. Note that the running of Yukawa couplings is also 
taken into account as well as the constraints coming from the top mass and possible b-r unification 
requirement. Finally, from this approach we get the allowed regions for the SUSY particle masses 
determined by the simultaneous effect of SUSY and GUT thresholds. From these results one 
obtains indications about the discovery potential of SUSY at LHG hence clarifying the role of 
future colliders. Note that the effect of possible intermediate thresholds could be also studied, but 
it would result in an extremely model dependent scenario. For this reason in order to be more 
predictive we prefer do not consider such possibility. 

The paper is organized as follows. In section II we give the generic details of SUSY-GUT model 
considered, whereas in section III we provide the Renormalization Group Equations for couplings, 
discussing the initial conditions and SUSY, GUT thresholds as well. Moreover this section contains 
a description of the numerical method adopted. In section IV we report our results whereas section 
V contains our conclusions. 


2 Overview on SUSY-GUT and naturalness principles 


Let us consider a generic SUSY-GUT based on a gauge symmetry G. Within the framework 
of TV = 1 supersymmetry, such a theory should contain vector (gauge) superfields V in adjoint 
representation of G, and some set of chiral superfields $ in different representations of G. A generic 
renormalizable Lagrangian can be written as^ 


>CsuSY = / + 


<feww+ / (few{<^) + h..c. 


( 2 . 1 ) 


where the first and second term yield the canonically normalized kinetic terms of the chiral su¬ 
perfields and the gauge superfields, and their gauge interactions. These terms do not involve any 

^For explicit notations, see e.g. Ref. [48]. 
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coupling constant apart of the gauge coupling, which is an 0(1) parameter. Third term describes 
mass and interaction terms between the fermionic and scalar components of chiral superfields $, via 
the superpotential hb(<f>) that is a general G invariant holomorphic combination $ containing the 
trilinear terms with 0(1) coupling constants and the bilinear (mass) terms. The soft supersymmetry 
breaking (SSB) terms can be presented in a form similar to Eq. (2.1) 


^SSB 


J d^ed^9p<l>h^^ + 


j d^epyvyv + 


d^9riW'{^) + h.c. 


( 2 . 2 ) 


making use of auxiliary superfields with non-zero F- and H-terms, respectively, rj = Mp9‘^ and 
p — M^9^9^, where the dimensional parameters Mp and Mp can be in principle different. In 
particular, Mp determines the size of gaugino mass terms, and via the third term in Eq. (2.2), 
contributes to the SSB terms for the chiral superfields. The D-term can be of the order of Mp 
once it is simply given by a direct product rjfj = Mp9^9^, while there can be also a direct H-term 
and in this case Mp ^ Mp. In general the function IE'($) should have the same structure of the 
superpotential 1T(<I>), but its couplings are not obliged to be the same of VE($). 

The chiral superfields can be divided into Higgs and fermion superfields distinguished by matter 
parity Z 2 , under which the fermion superfields change the sign while the Higgs superfields are 
invariant. In this way, the superpotential W has two terms related to the Higgs and Yukawa 
sectors, namely 

W = Wniggs + HVukawa • (2.3) 

The Higgs superpotential Hniggs is responsible for the v.e.v.’s breaking both the gauge symmetry 
G down to MSSM and then to SM. 


2.1 MSSM limit of a SUSY-GUT: mass scales and Higgs sector 

The MSSM represents the low energy limit of a generic SUSY-GUT. Here, we are going to discuss 
the different energy scales entering in the SUSY breaking pattern, paying particular attention to 
the Higgs sector. In MSSM the superpotential of Eq. (2.3) is explicitly given by 

IUmssm = , (2.4) 


which contains chiral superfields corresponding to three families of quarks and leptons (*, j = 1, 2, 3 
are family indices), and two Higgses Flu and Hd- The SSB F-terms repeat the structure of VEmssm, 
and contain also the soft Majorana masses of gauginos (bino, neutralinos and gluinos for a = 1, 2, 3) 


Cp = AlQ,u')Hu + A^QdlHd + Al^efl^Hd + pBpHuHd + . (2.5) 

Note that the dimensional parameters Bp and TOq are expected to be parametrically of 

the order oi Mp, though in the model dependent context the gaugino masses can be in principle 
different from other F-terms. On the other hand, the soft masses of all scalars including the Higgses 
are given by F-terms 


Cp) — j 


+ mlpufuj + miijdl 


2 JUJc 


■ wiy 


V. Li 


-mlfefe^^^+MlHlHu 


-MjWdHd. ( 2 . 6 ) 


All these dimensional parameters are of the same order of magnitude given by the energy scale Mp. 
These terms, in principle, can be parametrically larger than F-terms, though they can naturally 
be of the same order of magnitude. Therefore, all soft parameters can be divided in three classes: 
the soft Majorana gaugino masses 0(Mp), the soft masses of the scalars as squarks and sleptons 
0(Mp), and the so-called /r-term that determines the Higgsino masses and contributes to masses of 
scalar doublets Flu and Fid. GUT implies that all gauginos must have the same mass at the GUT 
scale (ttiq = me Va), and the similar mass unification can be assumed for masses of the squarks 
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and sleptons entering in the same GUT multiplet. Hence the differences between the true masses 
of squarks and sleptons, or between the masses of gluinos and neutralinos, are simply due to the 
running. In this way, SUSY-thresholds, namely the possible splittings among the masses of SUSY 
particles, can be defined by three mass scales only, which are gluino mass ifig 0{Mp), squark mass 
fhsq 0 {Md), and the Higgsino mass ttih 

A comment concerning the Higgs scalars and Hd is in turn. Their mass matrix 




V + ) 


(2.7) 


involves the mass parameters of three different origins: mass terms and Mj (SSB U-terms), 

(SSB U-term), and supersymmetric /u-term, which can be of different orders of magnitude. If 
the SSB parameters assume large values, namely larger than few TeV or even larger, one cannot 
pretend to recover the so-called natural solution, which would link the EW scale to the SSB scale. 
Hence a certain fine tuning condition has to be imposed. In particular, one eigenstate h should 
have a small negative squared mass 

-m^ = \ + Ml + Mj - ~ -(100 GeV)^ , (2.8) 

and that will be identified with the SM Higgs. It should get a v.e.v. {v = 250 GeV) breaking the EW 
symmetry and leaving behind the Higgs Boson with mass Ml = 2m? « (125 GeV)^. The other mass 
eigenstate h is heavy, with M| = 2/i^ -|- Ml + Ml + m^, and it should decouple at the SSB scale. In 
this way, the mixing angle between the Higgses Hu and Hd is given by tan2/3 = 2fiBg/{Ml — Ml), 
and hence the v.e.v. of h is placed between them Vu = v sin /3 and Vd = v cos /3. 

In principle, U-terms Mu^d could be larger than B^ and /r. However in this case the fine tuning 
of Eq. (2.8) would not be possible. On the other hand, a F-term B^ much larger than U-terms 
is not very natural. Moreover, for /r,i3^ Mu^d one would get tan ,5 « 1, which would make 
problematic to accommodate the 125 GeV Higgs. Therefore, to impose a realistic fine tuning only 
one possibility is left, namely to take soft parameters Mu^d^ and supersymmetric parameter /r 
all of the same order of magnitude, which looks rather embarrassing in the context of the generic 
SUSY models. However, in the framework of the SU{6) model [19] with pseudo-Goldstone solution 
for the gauge hierarchy and Doublet-Triplet splitting problems this situation arises rather naturally. 
In this case in fact, the supersymmetric /r-term emerges as a result of SUSY breaking, and at the 
leading order approximation the condition of Eq. (2.8) is straightforwardly obtained. 

In view of this, if one extends the needed conspiracy between the soft F- and Z3-terms in the 
Higgs sector to the squark and slepton masses (Z3-terms) and the gaugino masses (F-terms) as well, 
then all supersymmetric partners will be expected to have masses of the same order of magnitude. 

Let us remark also that for the proper condensation of the Higgs field a large quartic coupling 
constant ruling its self-interaction is needed, namely, one should have A = Ml/2v^ « 0.13. On 
the other hand, in the framework of the MSSM such a coupling is mainly provided by the gauge 
F-terms of SU{2) x U{\). Hence, at a scale of the order of the mass terms present in Eq. (2.7) 
one has A = ^{g^ + g' ) cos^ 2/3, where g and g' are the running gauge couplings of SU{2) and U{1) 
respectively. At the electroweak scale the gauge couplings are too small for saturating the needed 
value and one only gets A « 0.07. From this result it is derived the famous Higgs mass limit in the 
MSSM Ml = Ml cos^ 2/3 -I- radiative corrections, which makes difficult to recover the LHG Higgs 
mass Mh ~ 125 GeV and requires huge radiative corrections. However, the running of the Standard 
Model couplings shows that A fastly decreases with energy [49], and at a scale of the order of 10 
TeV it almost matches the needed value, unless tan/3 is too small. As it was shown in Ref. [50], by 
taking into account radiative corrections of reasonable size, the MSSM Higgs picture in the above 
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described decoupling limit is indeed compatible with the 125 GeV Standard light Higgs. This is the 
case for a scale of SUSY larger then few TeV, and if tan /3 is not too small, namely tan (3 > 2. 

Therefore, the Higgs physics, modulo specific conspiracies, indicates that the SUSY scale should 
be above few TeV, and similarly the bounds from the electric dipole moments settle the lower limit 
to few TeV as well. Generally flavor violations fix more stringent bounds, and for arbitrary soft 
masses and A parameters they would require SUSY scale larger than ~ 100 TeV. However, there can 
be (flavor) symmetry reasons that providing approximate alignment of the soft SUSY parameters 
with the Yukawa couplings reduce the lower limit to ~ 10 TeV. Moreover, in some theoretical 
scenarios based on flavor symmetries the flavor limits can allow SUSY scales as small as 1 TeV 
providing interesting relations between the soft mass matrices and SSB trilinear terms with the 
fermion Yukawa matrices as e.g. fhlij = m^dij + m\{Y^^Y^)ij + and A^j = msYfj, 

Wo, 1 , 2,3 being soft mass parameters [40-43]. Therefore, one can conclude that SUSY scale larger 
than few TeV is quite natural and if SUSY indeed lives near this lower bound, then it should be 
discovered in the LHG run H. 

2.2 A generic SUSY-GUT in the SU(5) bottleneck 

In minimal SUSY SU{5) [8, 9], the Higgs sector contains the following chiral superflelds: E in the 
adjoint representation 24 of SU{5) and H, H respectively in the fundamental representations 5, 5. 
In this case one has the following expression for the Higgs part of the superpotential in Eq. (2.3): 

lUHiggs = + + MhHH + ^ HYH . (2.9) 

There exists a supersymmetric minimum in which the adjoint held E gets a v.e.v. breaking SU{5) 
down to SU{3) x SU{2) x 17(1): 

(E) = Usxdiag(^^,^,^,-l,-l) , = (2-10) 

As a result, the “leptoquark" gauge bosons X,Y of SU{5) get a mass 

5 

Mg\jt = Mx,y = 2 \/27raGU7rUs , (2.11) 

where ogut is the SU (5) gauge coupling, while the Higgs supermultiplets in the broken phase have 
masses 

Me = -AeUe , Md = fJ- = Mh — Us , Mr = Mh + ~ 3^^^ ’ (2-12) 

In the previous expression Ms denotes the mass of the color octet and weak isospin triplet fragments 
of the SU{5) adjoint Higgs E ((8,1) © (1,3) under SU{^) x SU{2)), which are degenerate in mass. 
The quantity Mt denotes the mass of the Higgs triplet supermultiplets T, T contained in i7, H, and 
Md is the mass (/r-term in Eq. (2.4)) of their doublet fragments and Hd- The latter, which 
determines the Higgsino masses in the MSSM, should be small, while Mt should be large, order 
the GUT scale, in order to avoid a too fast proton decay mediated by the color triplet Higgs or 
Higgsino exchanges. In minimal SU{5) [8, 9], the price for such Doublet-Triplet (D-T) splitting is 
the fine tuning between the parameters in the superpotential of Eq. (2.9). In particular, Eq. (2.12) 
shows that two large values Mh, ^Ue 10^® GeV should be fine tuned for obtaining ^ 10^ GeV, 
with the precision of about 10“^^. However, this condition is sufficient when one takes into account 
the SSB terms in Eq. (2.2), in which the function rjW in this case can be presented as 

^W^Higgs = S' + ^ + Mh BhHH + ^ A^MEH^ , (2.13) 
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where A and B are dimensional parameters of the order of Mp. The presence of these terms 
shift the v.e.v. of S and generate its non-zero F-term, and at the end, the full account of these 
contributions generates the soft term fiB^ in Eq. (2.5), as fiB^ = — A^ — Bp + B^) + 

0{Mp). Therefore, for adjusting this value 0{VsMf) to the needed value ~ Mp ~ 1 TeV^, 
the parameter Ah — Ay: — Bh + By, which is 0(Mf), should be extremely fine tuned. Therefore, 
this situation cannot be called a natural solution of the hierarchy and D-T splitting problems. 

In minimal 51/(5) model the problem of D-T splitting has only a technical solution, given by 
fine tuning of fjL in Eq. (2.12) that is stable against radiative corrections. This situation is then 
worsened by the need of another fine tuning in the soft term ^B^. Of course, minimal 517(5) model 
is not realistic, and one should not be surprised to find that fine tunings are required. 

Below we describe more natural situations when fine tunings can be avoided in specific GUT 
models. Nevertheless, we shall consider minimal SU{5) as a prototype model, in the sense of SU{5) 
bottleneck condition, for understanding the possible threshold corrections.^ 

There are several realistic models in which the D-T splitting problem can be solved without fine 
tunings. In particular, in SU{5) this can be done via the “Missing Doublet Mechanism" (MDM) 
[10-14], in 50(10) via the “Missing v.e.v. Mechanism" (MVM) [15-18], while in SU{6) via the 
“pseudo-Goldstones instead of Fine Tuning" (GIFT) Mechanism [19-22, 25, 27] . In these models, 
the required patterns of the superpotentials are usually obtained by imposing some discrete symme¬ 
tries, which are guaranteed at the level of renormalizable couplings, but cannot provide suppression 
of dangerous high order operators destabilizing these solutions. It is interesting to note, however, 
that stable solutions at any order can be achieved by making use of the anomalous U{\)a gauge 
symmetry (see Refs. [14, 18, 24] respectively for the MDM SU{5), MVM 50(10) and GIFT 5C/(6)). 
In these models the soft parameters like By and Ay for the heavy GUT breaking superfields can 
be quite large, without creating additional fine-tuning problems. Let us briefly describe them here. 

In SU(5), the MDM [10, 11] contains the Higgs superfields in representations $ ^ 75, i7 ~ 5, 
H ~ 5, 4* ~ 50, ^ ~ 50, with the following superpotential terms: 

W = -t A4>^ -k Ml4-^ -k Aii7$^ -k A2i74>4' -k ^lHH , (2.14) 

with M and Mi being the mass parameters order Mqut and A’s being 0(1) coupling constants. 
SU{5) is broken down to 51/(3) x SU{2) x U{1) by the v.e.v. of $ ^ 75. The latter represents a less 
economic multiplet replacing the adjoint F of the minimal 51/(5), which also generates the mixing 
between the color triplet fragments in the Higgs 5- and 50-plets, whereas there are no doublets in 
the 50-plets. In this way, all color triplets are heavy, with mass of order Mgut, while the doublets 
in H,H remain light, with mass given by the /x-term. In this theory the new soft parameter 
(see Eq. (2.13) for its analogous definition) can be taken large without inducing a large B^. The 
unpleasant feature of this model is that the parameter /x should be taken very small with respect 
other masses M, Mi in an ad hoc way. 

Seemingly this model should satisfy our 5t/(5) bottleneck condition, since it is based just on 
51/(5). However, the situation is not so simple. The problem is that it involves the huge superfields 
75, 50 and 50, and if their fragments remain just slightly lighter than MguTj their large threshold 
corrections can completely ruin the gauge coupling unification. Moreover, the symmetry motivated 
versions of the MDM model [13, 14] contain twice as large amount of these fields. Therefore, in the 
context of these models the crossing of gauge couplings, with the todays precision, can be hardly 

^In particular, the masses My, 3-re important for our analysis since, if they are smaller than MquTi they 
would give significant threshold corrections near the GUT scale and affect the gauge coupling unification. However, 
within our paradigm of naturalness, couplings Ay; and g are assumed to be 0(1), and thus My, My > Mgut- 
Nevertheless, we study also situation when these fragments can be relatively light, taking a possibility that within 
the natural spread of 0(1) values, some of the above couplings could be as small as 0.1. 



considered stable for large threshold corrections. 


For SO(IO), in the MVM [15, 16], the philosophy is similar: the Higgs doublets remain massless 
since some of the GUT-breaking fields have vanishing v.e.v. along the direction that would give 
them a mass, whereas they couple with the triplets with non-zero v.e.v.. Therefore also in this case 
the protection of the doublet sector is due to group theoretical reasons. Therefore large soft terms 
^ 10 TeV in the heavy sector will not influence /i and B^. 

Regarding one step unification condition, it is not a priori guaranteed in 50(10). In fact this gauge 
group can be broken down to SU{3) x SU{2) x U{1) not passing through the SU{5) bottleneck, but 
e.g. via Pati-Salam subgroup 50(4) x 50(2) x 50(2) and subsequent symmetry breaking chains. 
This can be settled by ad hoc choice assuming that the Higgs superfields, which break 50(10) down 
to 50(5) subgroup, have the largest v.e.v.. However, there remains the problem that this model, 
and 50(10) models in general, are not very economic as far as the Higgs sector is concerned. The 
minimal set of GUT superfields necessary to achieve the correct symmetry breaking and MVM pat¬ 
tern within the renormalizable superpotential, contains one 54-plet, at least three 45-plets, 16 0 16 
and two 10-plets [17, 18[, whereas in some versions the much larger multiples as 210 and 1260126 
are used. Large threshold corrections due to these fragments can easily affect the crossing of gauge 
couplings and thus transform the gauge coupling unification from a prediction to an accidental fact. 

The SU{6) model [19-24] is based on very simple set of GUT superfields, in fact as simple as 
the minimal 5t/(5), SU(6) gauge symmetry is broken by two sets of superfields: two fundamental 
representations ~ 6 and H ^ Q that break SU{6) down to SU{5), and one in adjoint represen¬ 
tation E ~ 35, that leads to the breaking channel SU{6) —)■ 5U(4) x SU(2) x U{1). As a result 
the two channels together break the SU{6) gauge symmetry down to SU{3) x SU{2) x U{1). The 
simple assumption on which the GIFT mechanism is based is that the Higgs superpotential does 
not contain the mixed term HY,H, so that it has the form W = IU(E) + W{H,H), where 

1T(E) = W{H,H)=Y{HH-V^). (2.15) 

As a result the superpotential acquires an accidental global symmetry 5[/(6 )e x SU{6 )h, which 
independently transforms E, and H, H superfields. Then, in the limit of unbroken supersymmetry 
the MSSM Higgs doublets Hd appear as massless Goldstone superfields built up as a combina¬ 
tion of doublet fragments from E, and H, H, that remain uneaten by the gauge bosons. Therefore 
in this limit p vanishes exactly. 

Supersymmetry breaking terms like As, B-^ shift the v.e.v.’s and also give F-terms to them, 
therefore generating B^ term for the MSSM Higgses. However, since these terms also respect the 
global symmetry 5[/(6)s x SU{6)h, the mass matrix of the Higgses in Eq. (2.7) is degenerate, and 
thus one Higgs scalar (combination of the scalar components of and Hd) still remains massless. 
Thus, even with arbitrary that give p ~ ^ i?s, there is an automatic relation between p and 

B^ terms that guarantees that the determinant of Eq. (2.7) vanishes. This degeneracy is removed 
only by radiative corrections due to Yukawa terms that do not respect the global symmetry, and the 
resulting Higgs mass will be of the order of p and B^, given by the mismatch in their renormalization. 
Therefore, in the case of large Hs ^ 10 TeV we are still left with a ’’little” hierarchy problem of the 
electroweak scale stability against 10 TeV scale discussed in previous section. 

Due to economic Higgs sector and the way in which the GIFT mechanism works in 5U(6), the 
SU{5) bottleneck condition is satisfied in straightforward way. In a most elegant way this occurs in 
the context of model [24] with anomalous U{1)a symmetry where H,H v.e.v.’s emerge essentially 
at the string scale. Below such energy the theory becomes a minimal SU{5) with one adjoint 24 
and two fundamentals 5 and 5 plus one 5t/(5)-singlet, altogether composing a 35-plet of 51/(6) 
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that remains a global symmetry at this point. 


Until now we have not considered yet the fermion sector and the Yukawa part of the superpo¬ 
tential in Eq. (2.3). Coming back to a naive minimal SU{5) where fermions are allocated in the 
representations 10^ © 5^, i = 1, 2, 3, one can write the Yukawa couplings 


W'Yukawa = , 


(2.16) 


which after SU{5) breaking reduces to the Yukawa part of VEmssm reported in Eq. (2.4). However, 
one immediately encounters two problems: 

i) this predicts = Yfj and thus degeneracy of the Yukawa eigenvalues between down quarks 
and leptons of all generations at the GUT scale, yd,s,b = ye,iJ.,T- While for the third generation the 
yi^ = y^ unification works (almost) perfectly, for the lighter generations it is dramatically wrong; 

ii) in this way we must introduce small Yukawa couplings for lighter generations. In fact while only 
top quark has a large Yukawa j/t ~ 1 (for large tan/3, yi, and yr could be large as well), the rest 
of the Yukawa couplings must be <?; 1 in any case. This would contradict our assumption of no 
artificially small parameters in the mother GUT theory. 

Fortunately, both these problems can be solved at one shot in a rather natural manner, assuming 
that, by some symmetry reasons^, only the third generation gets masses from the renormalizable 
Yukawa terms of Eq. (2.16), while the masses of the first two generations are due to higher order 
operators 


H 

10,10, — + lOilO,-^ + a: 


' Mp 


^ Ml 


- EH 


6 ^ 10 , 5 , 


- E^H 


Ml 


+ .... 


(2.17) 


In the previous expression Mp denotes the Planck energy scale and i,j = 1, 2,3, with the coefficients 
a,b ^ 0{1) (see e.g. in [51]). The expression (2.17) after SU{5) breaking gives the Yukawa terms 
for the light families and their mixing with the third family. In this case we gain two things: the 
light fermion Yukawa couplings are no more degenerate, since they get Glebsch factors from the 
v.e.v. of S, and hierarchy between the families can be naturally understood in terms of small factor 
(E) /Mp ^ 10“^ ^ 10“^. The effective high order operators can be obtained by integrating out some 
heavy fermions in vector-like representations 5 © 5, 10 © 10, with masses 0{Mp), with which light 
chiral families 10i©5i mix with the 0(1) Yukawa couplings [52-55]. In this way, the mother theory 
becomes free of small couplings, namely all coupling constants are 0(1). It should be also noticed 
that due to these terms, the yt = yr unification in third family is not anymore exact but it gets 
corrections, typically yi, = ?/,-(l ± 0(y^/yr)) while in some predictive scenarios with asymmetric 
textures [18, 56, 57] corrections could be even larger. 

This situation emerges in a very natural way in the context of SU(6) model [19]. Its minimal 
anomaly free fermion content contains three generations of chiral superfields in multiplets that are 
a SU(6) decomposition of a 27-plet of Eq, 27 = 6 © 6 © 15, containing two (anti)fundamental 
representations, and a two-index antisymmetric one.® Their SU(5) decomposition reads 6 = 5 © 1 
and 15 = 5 © 10. The only possible Yukawa couplings are gikl^fikH, i = 1, 2, 3 and fc = 1, 2, ... 6 , 
where the matrix gn- can be taken block diagonal without lose of generality, e.g. with non-vanishing 
elements being only gi 4 , 525 ) 536 - After symmetry breaking SU(6) —>■ SU(5) by the v.e.v. of H 
(Vh), these couplings combine three fragments 5^ C 15^ with three combinations of six 5 fragments 


^Typically due to family symmetries (in SU{6) model this occurs automatically due to pseudo-Goldstone nature 
of the Higgs superfields) only top quark can get the mass from renomalizable coupling with yt ~ 1, while other 
masses come from higher order terms suppressed by Mp. It is worth to notice, however, that in spite masses of b and 
T emerge from higher order operator, j/j = j/r still holds with the modulo corrections related to hierarchy between 
the SU{6) 517(5) and 517(5) -)■ 517(3) x 51/(2) x 1/(1) breaking scales with Vs/Ve ~ lO"! -b [22, 23]. 

®For a comparison, 50(10) decomposition of the Eg 27-plet is 16 © 10 © 1 containing the spinor representation 
16 , but also additional non chiral fragments 10 © 1 . 
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from 6fe into massive particles with M ~ Vh- Consequently, theory reduces to SU{5) with three 
generations 10^ © 5^, i = 1,2,3. However, non-renormalizable Yukawa couplings can be written 
that could generate fermion masses. 

It is possible to make a minimal modification of the model by introducing a superfield 20 
(three-index antisymmetric representation) [22, 23] whose SU{5) decomposition is 20 = 10 © 10. 
Since 20 is a pseudo-real representation of SU{6), its mass term is not allowed: convolution of 
20^ with a Levi-Chivita tensor is vanishing. Therefore, its 10 © 10 fragments remain massless 
before SU{6) symmetry breaking. In this case, the Yukawa couplings A 20 20 E and A'2015^1/ 
that explicitly violate the global SU{6)^ x SU{6 )h symmetry can be introduced. The latter one 
combining 10 C 20 with 10 C 15 provides massive states, while the former one reduces to a SU{5) 
Yukawa coupling for the fragment 10 C 20 with 5 fragment of E, A 10 10 5e. In this way, an upper 
quark from 10, to be identified with top quark, gets Yukawa coupling with the pseudo-Goldstone 
fragment Hu C 5e, with yt ^ A ^ 1. The t-quark is the only particle getting in this way mass 
0(100) GeV [22, 23[. Masses of other particles should necessarily emerge from the higher order 
operators. Namely, masses of b and r come from the operator (Xb/Mp) 20 implying yi, = yr 

modulo small corrections, namely j/b = j/r(l + 0{V^/Vh)') ■ The mass of c-quark is induced by 
the term (Ac/Mp) 152l52i?^E, with a nice implication j/c 2/6,r ~ (I4r/Mp)^ pointing towards 
small or moderate tan/?. It is worth to remark that this feature also implies the lower bound 
(Vh/MpY > 0.01, i.e. Vh Y 10^^ GeV. Hence, Vh ^ Vs and the SU{b) bottleneck condition is 
naturally satisfied. 

In this framework, the SSB D-terms of Eq. (2.16), which generate masses for squarks and 
sleptons, are given by 

Z/Sihs) ^ ^ ^2^^ gtg^., (2,18) 

In our analysis, we allow the mass matrices mio and 7715 to be different. Different squark and 
slepton fragments, which have the same soft masses in the GUT limit, split at low energies due to 
the running. The effect of such difference on the compatibility of the model has been considered 
in the following analysis. Note that fh^o and ms would be exactly equal in S'O(IO) inspired 
models at Mqut- As for inter-family splitting of soft masses, as well as for healthy pattern of the 
trilinear terms in Eq. (2.5), they can be alligned with the Yukawa terms through MFV relations 
between the Yukawa couplings and soft parameters for squarks and sleptons, as e.g. m-Aj = 
rriQSij + m^(Y“fY“)y + m2(V"fY“)f^ and H"- = mVVpj, where moq,2,3 are soft mass parameters. 
Therefore, due to small values of the Yukawa couplings, mass spectrum of the squarks and sleptons 
must have inter-family degeneracy with possible exclusion of stop since yt 1. In any case, this 
unification of fermion and sfermion masses and interaction patterns provides a chance for a natural 
suppression of SUSY induced flavor violating effects, and allows to lower the SUSY scale down 
to 1 TeV. It is important that such relations can be obtained also in the context of realistic and 
predictive SUSY-GUT scenarios for the fermion masses and mixing [40-44] . 

3 Renormalization Group Equations: initial conditions and threshold 
corrections 

Let us consider the Renormalization Group Equations (RGEs), up to 2-loop order, for a set of 
couplings Xi that can be generally cast in the following form 

= + (3.1) 

where t = ln{M/Mo), and Mq represents the renormalization energy scale at which we impose the 
initial conditions. The set of running couplings and the corresponding expressions for and /?^.^ 
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Figure 1. Running of the gauge couplings (blue), (green) and (red) with Msusy = 2 TeV up 
to 2-loop in /3-functions. For each a~^ it is reported the 3a error band due to the experimental uncertainty 
at Mz- The range for Msusy compatible with the measurements within Scr is from ~ 1 TeV up to ~ 3.5 
TeV. The inset is a focus of the crossing region. 

depend on the energy regime considered, simply denoted as SM (non SUSY) or MSSM (SUSY), and 
on the massless particle content at that particular energy scale. The complete set of SM RGBs is 
provided by [58-63] , while the MSSM RGBs are reported in [64] . We properly assume that only the 
Yukawa couplings for the heaviest particles of each family generation (i.e. top and bottom quarks 
and tau lepton) give a considerable contribution. Moreover, in the running of gauge couplings we 
consider only the top Yukawa contribution since the others are negligible at all energy scales for 
moderate tan /3. Moreover, in the MSSM /3-functions we neglect the terms related to trilinear soft 
couplings of Bq. (2.5). 

As usual, starting from low energy, it is convenient to fix the renormalization mass Mq at the 
BW scale Mz = 91.1876 ± 0.0021 GeV [34], where the gauge couplings at = gfl^iT result to be in 
the MS renormalization scheme [34] 

a-i(Mz) = 59.008±0.008, a^^(Mz) = 29.569± 0.014, a3(Mz) = 0.1184± 0.0007, (3.2) 

where f = 1, 2, 3 stands for U(l), SU{2), SU{3), respectively. As it is well known, in the running 
from low to high energy scale we have to add every particle above its individual threshold given by its 
mass, and to impose at the same time corresponding matching conditions. In the SM this happens, 
for instance, in the case of Higgs (Mh = 125.7±0.4 GeV) and of top quark {Mt = 173.21±0.51±0.71 
GeV) [34] . 

Starting from Mz, with the increase of energy we have at a certain scale the transition from 
SM to MSSM. Here, we have two possibilities, namely, the single-scale and multi-scale approaches 
[65] . In the single-scale approach, the transition between SM and MSSM occurs at a given effective 
energy scale simply denoted as Msusy- In this scenario all SUSY particles share the same mass. 
In Big.l one shows the gauge couplings unification in a SUSY scheme with a single effective energy 
scale Msusy ~ TeV. In achieving such result, as it is well known, 2-loop /3-functions play a relevant 
role in increasing the value for A/susy needed. However, the simplifying assumption of single-scale 
approach is very unnatural, since, in principle, the SUSY particles mass spectrum might be not 
degenerate. When such occurrence is considered we have the so-called multi-scale approach that is 
adopted in our analysis. 

In the framework of SUSY-GUTs under the SU(5) bottleneck ansatz, sparticles, even though 
with the same mass at Mgut, naturally acquire different masses al low scale due to the renormal¬ 
ization group equations. In our model, according to Bqs. (2.5) and (2.18), at Mqut there exist 
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three parameters related to SUSY soft masses tog (SSB F-terms) for gauginos and, toiq and ms 
(SSB D-terms) for matter sparticles. Let us denote with fhg, tow, hT-sq and TOsi the pole masses 
of gluinos, neutralinos, squarks and sleptons at low energy scale, respectively. Due to the SU(5) 
bottleneck ansatz, at Mgut we can assume TOg/mw = 1 and rnsq/ifisi = miotni^ = 1. It is worth 
while noticing that the simplicity assumption rhiojrhs = 1 does not sensibly affect the final results 
since the main contribution to gauge coupling running comes from gauge and Higgs sectors as will 
be clarified in the following. Finally, there exists also another parameter Wh that is the mass of 
two higgsinos of the same order of supersymmetric parameter /i. Note that the mass of the second 
Higgs is of the same order of the masses of matter sparticles. 

Let us define and Mmax as the minimum and the maximum of SUSY particles mass 

spectrum, respectively, in a given model. Concerning the SUSY thresholds we adopt the following 
simplicity ansatz. We use below Mmin 2-loop SM RGBs, while above Mmax 2-loop MSSM RGBs. 
Within these two scales, we apply 1-loop MSSM RGBs for each new SUSY particle, i.e. adding the 
contributions of each particles at the corresponding mass (see Ref. [66]), while we use 2-loop SM 
RGBs for SM particles. 

Moreover, since the SM /3-functions are typically evaluated in the MS renormalization scheme, 
while the MSSM ones are obtained using the DR, one should take into account matching relations 
between the couplings corresponding to the different schemes. However, since the transition between 
MS and DR gives a shift in the unknown SUSY particles masses smaller than 1%, for the level of 
accuracy of the present analysis we can safely neglect such effect. 

For Yukawa couplings of the corresponding particles, we use 2-loop SM RGBs up to Mmax and 
2-loop MSSM RGBs above such scale, imposing the suitable matching conditions 

^5M^^MSSMsin/3, COS fi, y^^ = y^^^^ COS /3 . (3.3) 

Finally, GUT particles which are lighter than leptoquarks {X, Y) modify RGBs before gauge 
couplings unify at Mgut- For instance, the mass of the (8,1) 0 (1,3) fragments of SU{5) adjoint 
Higgs, Ms, or of triplet fragments (3,1) 0 (3,1) in H,H, could be below Mgut if the coupling 
constants As or ^ in Bq. (2.9) are small. Indeed, by comparing Bqs. (2.II) and (2.12) we can define 
the following parameters 

Mgut V27raGUT Mgut v^Sttogut /„ .x 

Xs = — ^ = - r -, XT = = -T-, (3.4) 

Ms As Mt s 

If the parameter ys is large, then the GUT threshold corrections due to fragments of B will be 
relevant. According to our naturalness conditions, we assume that As is not unnaturally small, 
and impose the constraint ys < 10. The same can be applied also for color triplet fragments T, T. 
However, the proton stability against dimension-5 operators require Mt > Mgut- Therefore, ^ ~ 1 
is favored, and threshold corrections related to xt are expected to be irrelevant. It is interesting 
to notice that in the context of SU{6) model [19], one automatically gets Ms = Mt as far as 
they emerge from the same adjoint Higgs E ^ 35. Hence, the requirement of the proton stability 
against Higgsino mediated dimension-5 operators, Mt > Mgut> also implies that there should be 
no threshold corrections neither from other fragments the superfield E. 

The b-T unification at GUT energy scale is also studied, by imposing at Mgut 

J/b(MGUT) = //.(Mgut) (l + O • (3-5) 

Note that since the ratio between two Yukawa couplings is essentially unchanged by the running, 
we have used their values measured at BW scale in order to estimate O [y^/yr)- The masses are 
Wb(TOb) = 4.18 ± 0.03 GeV (evaluated in MS scheme), Mr = 1776.82 ± 0.16 MeV and M^ = 
105.6583715 ± 0.0000035 MeV [34] . 
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Figure 2. Steps of the iterative method implemented in the Mathematica program which solves the RGEs, 
imposing gauge unihcation at Mqut- All red solid lines mean that the differential equations are resolved 
from Mz up to Mqut, whereas the black dashed line means that the corresponding RGE is resolved from 
high to low energy scales. The circles represent the initial conditions imposed at Mz, while the black 
triangles denote the output parameters that must be compatible with experimental bounds. The energy 
scales, Mmin and M^ax, are respectively the minimum and the maximum of SUSY particles masses in a 
given model. The quantity Ms denotes the mass of E Higgs fragments (8,1) © (1, 3). Note that all these 
energy scales dehne the SUSY and GUT thresholds. 


We have developed a Mathematica program® which resolves all the RGEs with numerical itera¬ 
tive method and takes into account all the matching and threshold relations discussed in the previous 
sections. The input parameters of the program are {fhh, rUg, TOgq, xs, tan,5} (hereafter simply de¬ 
noted as a model), whereas the outputs are {q:3(Mz), Mqut, ogut, t/t(MGUT), yb(MGUT), yT(MGUT)}. 
The quantity xt has not been considered in the analysis since it contributes in an almost irrelevant 
way ^ 1) to the running. 

Since the output quantity 03(Mz) must be compatible with its experimental value it represents 
a compatibility constraint that has to be satisfied by the particular model chosen. In other words, 
scanning on the possible input values one has to discard those choices producing a value of 03(Mz) 
outside the experimental range. This yields, at the end of the procedure, a set of compatible 
SUSY-GUT models. 

Another compatibility constraint concerns Mqut since it is straightforwardly related to the 
proton lifetime Tp. In a generic SUSY-GUT scheme, the proton decay is induced by dimension-6 
operators (mediated by gauge bosons X,Y), and by any type of effective operators allowed by 
a general operator analysis [70, 71]. In this study we focus our attention on dimension-6 oper¬ 
ators responsible for the relevant decay channel p —>■ e^7r®, namely the exclusive decay channel 
posing the most stringent bounds. The corresponding experimental bound is T”^ (p —>■ = 

Tp/Br (p —7> e+TT®) > 1.29 • 10®"^ yr [72] with 219.7 kt-yr of data at 90% confidence level. From Ref. 
[73] one gets for the partial width 


r(p^e+^®) - + + [(^r)) + (^r)) (^ + 


. (3.6) 


By substituting in the above expression: = 0.13 GeV, an = —0.0112 GeV®, D = 0.8, F = 0.47, 

« 2.5 and « 2.6 (see Ref. [73] for the definition of the parameters) and using [34] one gets 


E"^ (p-^ e+TT®) « 2.4-10' 


32 f Mqut 


1016 GeV y 


^ yr > 1.29-10®'‘yr 


Mgut >3.ioi6GeV. (3.7) 


U«GUT 


60ther programs able for such a purpose are for instance RunDec [67], SPheno [68] and SoftSusy [69] . 
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In the following we adopt this conservative limit that does not depend strongly on the details of the 
fermion masses and mixing. It is interesting to remark that the new effective operators of dimension- 
4 and -5 can be induced by the SSB terms when the heavy gauge bosons of GUT are integrated 
out [74, 75]. After being dressed by gauginos, they transform into dimension-6 operators, and in 
the case of large soft parameters, as it can be in our case, the proton lifetime in gauge mediated 
channels can be strongly affected (enhanced or even suppressed) depending on the pattern of soft 
terms in the heavy Higgs sector [75]. 

The dimension-6 operators induced by the exchange of heavy color triplet scalars in i7, H are 
suppressed by small Yukawa couplings of the first generation fermions, even though in some models 
with flavor symmetry [76, 77] their contribution can be dominant. In generic SUSY-GUTs, the 
dominant contribution to proton decay usually comes from dimension-5 operators induced by the 
color triplet Higgsino exchange [78, 79] . After being dressed by gauginos, they induce proton decay 
dominantly via the channel p —>■ K^V. The experimental limit on the latter is quite stringent, 
namely rp/Br(p— K'^V) > 5.9 • 10^^ yr at 90% confidence level [80] . This bound excludes the 
naive supersymmetric SU(5) model [81, 82], which is however already excluded by the wrong 
prediction for the light quark masses. For models reproducing realistic fermion masses the impact 
of these dimension-5 operators is strongly model dependent. In literature have been discussed many 
examples of models in which dimension-5 operators can be suppressed by particular symmetry 
reasons [13, 83-87]. By taking the constant ^ ~ 1, the color triplet Higgsinos move to mass scales 
higher than the GUT scale. In addition, having squark mass large enough, say 10 TeV, the Higgsino 
mediated dimension-5 operators indeed become safe. 

As discussed before, the running of all couplings is evaluated using 2-loop /3-functions, while the 
threshold effects are studied at 1-loop level only. Since the 2-loop RGBs are an extremely involved 
set of coupled equations, we employ an iterative method whose steps are shown in Fig. 2. In order to 
impose the gauge unification at Mgut, for a given choice of input parameters and starting from the 
known values of Mz, Mh, Mt, Mb, Mt, a)"^(Mz) and previously reported, we iteratively 

• evolve ai and a 2 , from Mz to higher energy scales; 

• find their intersection point that defines the values of ogut and its corresponding unification 
scale Mgut; 

• evolve backward as using (aGUT,M gut) as initial point, obtaining as{Mz)'i 

• starting from this solution for running gauge coupling constants we introduce the contribu¬ 
tions of ?/t) Vh and Ut in such a way that they consistently reproduce the known values of 
Mt, Mb and Mr- This fixes at the end the corresponding values for j/t(MGUT), yb(MGUT) 
and yr(MGUT)- Note that the other Yukawa couplings are not considered since they give 
subdominant contributions. 

This procedure must be iterated until it converges. The convergence conditions is that the absolute 
difference |a3^^(Mz) — a\{Mz)\ < 10“®. Such precision is sufficient since it has to face an experi¬ 
mental uncertainty of the order of 10“^. We have also verified that the modeis considered remain 
perturbative during the running of all couplings. 

In the multi-scale approach, there exist in principle five masses parameters, i.e. mg, mwj di-sq, 
rhs\ and mh, which define the SUSY thresholds. In order to estimate the values for the two masses, 
mw and msi, assuming rug/mw = 1 and rhsq/rhs\ = rhxo/rhs = 1 at Mgut, we resolve backwards 
the RGBs of gauginos and sparticles soft masses in the single-scale approach, where Msusy is of 
the same order of magnitude of rhg and mgq, respectively. In particular, concerning matter SUSY 
particles we assume rngq/rngi = mq/mb at MgusY = disq- 

Goncerning the GUT threshold, for a given value of xs, we iteratively find the value of Ms 
through which the chosen value of xs is achieved. 
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Parameters 

Ois,{Mz) 

Mgut 

xe: 1 ^ 10 

+ 2.0 % 

+106 % 

rhg: 1 —>• 10 TeV 

- 2.2 % 

-47 % 

rhsq'- 1 —>■ 10 TeV 

+ 0.2 % 

- 2.6 % 

rrih: 1 —)■ 10 TeV 

-4.5 % 

-16 % 


Table 1. Level of dependence of a^^Mz) and Mqut on the main parameters related to SUSY and GUT 
thresholds. The dependence on tan /3, at this level, results to be negligible. 

4 Results 

In our analysis we try to determine, among the compatible SUSY-GUT models, namely the possible 
choices {thh, mg, rfisq, xs, tan/3}, the energy scale above which SUSY signatures have to show up. 
To define such a scale, let us consider a given compatible model where the corresponding SUSY par¬ 
ticles have masses mn, m-g, and mgq that admit a minimum. By scanning on all compatible models 
for fixed X's we can determine the maximum of previous minima (once that we have marginalized 
with respect to tan/3) hereafter denoted by 3 Wub(xe)- 

In Table 1 we report the behaviors of 03 (Mz) and Mqut as a function of the SUSY and 
GUT thresholds. In particular, concerning 03 (Mz) one observes two opposite dependences on the 
thresholds. The increasing in the GUT threshold, X's, produces an analogous increasing in a 3 {Mz), 
whereas the opposite occurs for SUSY thresholds. This can lead to a possible balancing between 
these different effects that allows for larger masses in the SUSY spectrum. As it appears from 
the Table, the squarks provides a negligible contribution to 03 (Mz) hence justifying the simplicity 
ansatz we already quoted mio/ms = 1 at Mgut- The behavior of Mgut on the thresholds is 
qualitatively similar, but the strong dependence on xn makes more difficult a possible balancing. 

In Fig. 3 we report on the first row of the panel the compatibility region for the masses mh, 
rhg, rhsq for different values of GUT threshold As it is clear, 771 ^ and rhg result to have a 
reasonable inverse correlation and this occurs almost independently of rhsq. From the Figure one 
can observe the presence of an upper bound for rhh that is smaller than the possible ones for the 
other particles, and that such value appears to be a monotonic increasing function of xs- Note that 
for simplicity in these pictures we do not superimpose the constraints coming from proton lifetime 
and b-T unification. 

The other rows of Fig. 3 provide the allowed regions of the input parameters once one has 
reduced the level of arbitrariness among the values of SUSY masses. In particular in the second, 
third and fourth row one assumes the relations rhh = higi, rhh = fh-yv and rhg = rhgq, respectively. 
Such situations cannot be straightforwardly derived by the 3d plots of the first row, and hence they 
complement such information. Moreover, these three different situations are in agreement with the 
fine tuning of Eq. (2.8), which require the masses rhg (SSB F-terms), rhgq (SSB D-terms) and rhh 
(supersymmetric ^-term) to be of the same order of magnitude. The superimposed dashed lines 
correspond to given values of Mgut, expressed in 10^® GeV, whereas the red line, if present, bounds 
from above the region allowed by proton lifetime. Note that Eq. (3.7) bounds the compatibility 
regions for small values of xe only. This is due to the fact that larger values of xs are able to 
spread and shift up to larger values the compatibility regions for SUSY particle masses (see second 
and third column of Fig. 3). Moreover xe affects the Mgut allowed region as well, as it is clearly 
shown in Fig. 4, where it can be seen how increasing the values of xe one gets larger values for 
the unification scale. The red-dashed regions of Fig. 4 are excluded by proton lifetime constraint 
ofEq. (3.7). 

The anti-correlation showed by some of the mass parameters suggests that to look for Mub(xe), 
it is convenient to assume particular relations among SUSY particle masses. In particular, in the 
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Figure 3. In the first row of the panel the compatibility regions for the masses mt, rhg, and rhaq for 
different values of xs are reported. The second, third and fourth row provide the allowed regions once 
they have been assumed the relations mj, = ttIsi, rhh = mw and thg = rhaq, respectively (green and yellow 
regions correspond to 99% and 95% CL respectively). The dashed black lines correspond to given values 
of Mgut expressed in 10^® GeV, whereas the red line, if present, bounds from above the region allowed by 
proton lifetime (Eq. (3.7)). 
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Figure 4. Allowed regions in the Mgut-ugut plane obtained for different values of GUT threshold xs- 
The red-dashed regions are excluded by proton lifetime constraint of Eq. (3.7). 






(a) (b) 


Figure 5. (a) Allowed region in the plane tan/3-i/t(AfGUT). The yellow (green) region corresponds to 

compatibility within 2a {3a) uncertainty. The blue, purple and red lines bound from below the allowed 
region once that exact or partial b-r unification at Mgut is respectively assumed (see Eq. (3.5)). (b) 

The quantity Mub as a function of xe. On the upper part of the picture, we report the values of As 
corresponding to the values of xe according to Eq. (3.4). 


TOh-TOg plane one can see that if a mass is very light, the other one has to be very heavy, whereas 
the quantity rhsq is not bound at all. This suggests that Mub(xi:) can be found by imposing 
TOh = Wg = rhsq in the allowed region, and looking for the maximum of such values. 

In the left panel of Fig. 5 we provide the allowed region in the plane tan j3-yt{MG\jT)■ As it is 
clear from the plot, for moderate tan/3 {naturalness requirement, namely small fine tuning in the 
Higgs and Yukawa sector) ytiMovT) is bound to be larger than 0.5. In the case of exact {yi, = y-r) or 
partial (j/f, = 0.90 yr and yb = 0.85 b-T unification at MguTj according to Eq. (3.5), the quantity 
//((Mgut) cannot be lower than ^ 1.5, ^ 0.7 and ~ 0.5, respectively. It is worth observing, that the 
b-T unification at GUT energy also imposes tan/3 < 3.0 for the case yb = 0.90However, larger 
values for tan/3 can be obtained by relaxing the relation of Eq. (3.5). In particular, the allowed 
region in the plane tan,d-yt (Mgut) is not bounded at all for yb = 0.80y^, which is considered as 
the case “without b-r unification". 
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In the right panel of Fig. 5 we show the quantity Mub(xs) as function of GUT threshold xy.- 
On the upper part of the picture we report the values of As corresponding to the particular values 
of xs- As it can be easily seen, too large values of xt. would correspond to excessively small values 
of As that can be considered unnatural. Moreover, to assume large values for xt. > 10 would also 
imply Mgut unnaturally approaching Planck scale. For these reasons in our analysis we bound 
the values of xe in the conservative interval 1 10. As shown in the plot, 6-r unification at Mqut 

prefers lighter SUSY particles, implying lower values for Mub(xe)- 

For all these reasons, one gets that the upper bound for Mub(xe) is ~ 20 TeV whenever 
Xs < 10. For energy scale larger than such upper bound a generic SUSY-GUT model satisfying our 
requirements has to show up via the detection of Higgsinos or gluinos. Such a prediction cannot be 
extended to squarks that could be much heavier. 

Note that, in case we had used more stringent bounds due to dimension-5 operators mediating 
proton decay (even though highly model dependent), we would have excluded larger portions of the 
allowed regions reported in Fig. 4 thus favoring larger values of xs- In this case, the determination 
of the upper bound of ~ 20 TeV for Mub would remain completely unchanged since it only depends 
on the largest values of xs compatible with our requirements. 

5 Conclusions 

Besides the recent detection of a new particle compatible with the Higgs boson that completes the 
spectacular set of experimental evidences supporting the SM, it is hardly to believe that a QFT based 
on the gauge group SU{3) x SU{2) x 17(1) can represent the deepest description of fundamental 
interactions. This consideration mainly follows from the observation that several phenomena or 
open problems suggest the presence of physics beyond the SM. Among them, the hint of unification 
of all gauge couplings for extreme large energy, and the problem related to a natural separation of 
very different energy scales in a field theory with scalars {hierarchy problem) strongly indicate the 
need for more profound schemes. 

The GUT paradigm once implemented in a SUSY framework is able to simply address all these 
problems. Nevertheless, since the 8 TeV LHG run I has greatly constrained SUSY-GUT, one can ask 
if such simple scenario is still viable. In the present analysis, by considering a generic SUSY-GUT 
obeying quite general assumptions {SU(5) bottleneck, consistency with third family fermion masses 
and with the experimental limit on proton decay; absence of special fine tunings among couplings, 
which are 0{1) at Mqut) we have analyzed the limits on the mass spectrum of SUSY particles 
once the experimental constraints at Mz and proton decay limit are considered. 

Parametrizing the SUSY mass spectrum in terms of three quantities related to SUSY soft masses 
rhg (SSB F-terms) for gluinos, rhsq (SSB D-terms) for squarks, and Wh (supersymmetric /u-term) 
for Higgsinos, we have looked at those values compatible with the previous requirements. The aim 
was the determination of the upper bound for the minima of compatible spectra. Such energy scale 
defines a threshold above which SUSY signatures have to show up. This study provides indications 
about the chances to unveil new physics on LHG run II and on future colliders. 

In order to be more conservative in our prediction we have discussed the general possibility 
that SUSY-GUT posses several thresholds. This means to assume possible splittings between the 
masses of SUSY particles (multi-scale approach), and analogously, to admit that some GUT particles 
can take mass below the unification scale (GUT thresholds). In comparison with the single-scale 
approach where the SUSY breaking scale is of order of few TeV, in the multi-scale scenario we find 
that SUSY particles masses can be much heavier. 

We claim that if a SUSY-GUT model is the proper way to describe physics beyond the SM 
and under the natural assumptions reported above, the lightest gluino or Higgsino cannot have a 
mass larger than ~ 20 TeV. The requirement of b-r unification at GUT energy scale slightly reduces 
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such upper bound of few TeV. It is worth observing that our conclusions are strongly affected by 
the experimental uncertainty on a 3 {Mz) and on the proton lifetime lower limit. According to these 
results a new generation of colliders able to achieve the 100 TeV [46, 47] energy range would have 
the chance to cover almost all the parameter space beneath a natural unifying scheme. 
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